function ims=drawbox(im,cen)
% usage: img = drawbox(img, cen_prev_mat);

if size(im,3) == 1 
    ims=reshape([im,im,im],[size(im,1),size(im,2),3]);
else
    ims = im;
end

cen = round(cen);

for i=1:size(cen,1)
    if isstruct(cen)
        [suby,subx]=ind2sub(size(im),cen(i).pidx);
        xmin=min(subx); xmax=max(subx); ymin=min(suby); ymax=max(suby);
    else
        xmin=cen(i,1); xmax=cen(i,1);
        ymin=cen(i,2); ymax=cen(i,2);
    end
    
    if xmin-10<=0 || xmax+10>size(im,2) || ymin-10<=0 || ymax+10>size(im,2)
        continue;
    end

    ims(ymin-10:ymin-6,xmin-8:xmax+8,1)=180;
    ims(ymax+6:ymax+10,xmin-8:xmax+8,1)=180;
    ims(ymin-8:ymax+8,xmin-10:xmin-6,1)=180;
    ims(ymin-8:ymax+8,xmax+6:xmax+10,1)=180;
    
    ims(ymin-10:ymin-6,xmin-8:xmax+8,2)=200;
    ims(ymax+6:ymax+10,xmin-8:xmax+8,2)=200;
    ims(ymin-8:ymax+8,xmin-10:xmin-6,2)=200;
    ims(ymin-8:ymax+8,xmax+6:xmax+10,2)=200;
    
%     if xmin-5<=0 || xmax+5>size(im,2) || ymin-5<=0 || ymax+5>size(im,2)
%         continue;
%     end
%     ims(ymin-5:ymin-5,xmin-5:xmax+5,1)=245;
%     ims(ymax+5:ymax+5,xmin-5:xmax+5,1)=245;
%     ims(ymin-5:ymax+5,xmin-5:xmin-5,1)=245;
%     ims(ymin-5:ymax+5,xmax+5:xmax+5,1)=245;
    
end
end